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Numerical demonstration of fluctuation dynamo at low magnetic Prandtl numbers 
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Direct numerical simulations of incompressible nonhelical randomly forced MHD turbulence are 
used to demonstrate for the first time that the fluctuation dynamo exists in the limit of large 
magnetic Reynolds number Rm S> 1 and small magnetic Prandtl number Pm <C 1. The dependence 
of the critical Rm c for dynamo on the hydrodynamic Reynolds number Re is obtained for 1 < Re < 
6700. In the limit Pm <C 1, Rm c is about three times larger than for the previously well established 
dynamo at large and moderate Prandtl numbers: Rm c < 200 for Re > 6000 compared to Rm c ~ 60 
for Pm > 1. Is is not as yet possible to determine numerically whether the growth rate of the 
magnetic energy is oc Rm 1 ' 2 in the limit Rm — > oo, as it should if the dynamo is driven by the 
inertial-range motions at the resistive scale. 

PACS numbers: 91.25.Cw, 47.65.-d, 95.30.Qd, 96.60.Hv 
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Introduction. The amplification of magnetic field by 
turbulent fluid motion, or dynamo, is believed to be the 
cause of cosmic magnetism [H, 0, Q . Two types of turbu- 
lent dynamo should be distinguished. The first is the 
mean-field dynamo defined as the growth of magnetic 
field at scales larger than the outer (energy-containing) 
scale L of the turbulent fluid motion. The second, which 
is the focus of this Letter, is the fluctuation dynamo (or 
small-scale dynamo) defined as the growth of magnetic- 
fluctuation energy at or below the outer scale [25j |. 

Reflcxively, one tends to think of turbulence as an 
effective mixing mechanism rather than a constructive 
agent. It is then a remarkable idea that random motions 
can amplify a magnetic field. The currently accepted 
qualitative explanation of how the fluctuation dynamo is 
possible is based on the notion of the random Lagrangian 
stretching of the field by the fluid motion 0, H, 
This picture depends on the assumption that the scale 
of the stretching (which is the viscous scale l u , because 
the motions there have the largest turnover rate) is larger 
than the scale of the field that is stretched (the resistive 
scale l n ). Whether this is true depends on the magnetic 
Prandtl number, Pm = Rm/Re, which, in most natural 
systems, is either very large or very small. When Pm 3> 
1, we have indeed l v /l„ ~ Pm -1 ' 2 <C 1 (see 0| and ref- 
erences therein). In contrast, when Pm <g 1 (while both 
Re > 1 and Rm > 1), one expects l v /l u ~ Pm~ 3/4 > 1 
Q, so the resistive scale is in the middle of the inertial 
range, asymptotically far away both from the viscous and 
outer scales. Can the field still grow? In the absence of a 
mechanistic model of the field amplification, the problem 
is quantitative: the stretching and turbulent diffusion are 
of the same order at each scale in the inertial range, so 
which of them wins is not obvious. For over 50 years, the 
resolution of this problem has fascinated and confounded 
several generations of scholars, who were intermittently 



convinced that the fluctuation dynamo did or did not 
exist. 

With the advent of modern scientific computing, nu- 
merical simulations have been used to build a case for 
dynamo as a very generic property of random and chaotic 
flows @, M, E3, EH, El- This case has recently been 
strengthened by a successful laboratory dynamo in a ge- 
ometrically unconstrained turbulence of liquid sodium 
[l3j ]. However, both in the computer and in the lab- 
oratory, it is nearly impossible to access the values of 
Re and Rm that are sufficiently large to resemble real 
astrophysical situations. It has been especially diffi- 
cult to model the case of Re ^> Rm ^> 1, correspond- 
ing to the limit of low Pm. The computational chal- 
lenge in this limit is to resolve two scale separations: 
L 3> In 3> l v . The low-Pm limit is encountered in the 
liquid-metal cores of planets (Pm ~ 10 -5 ), the solar con- 
vective zone (Pm ~ 10~ 7 — 1CP 4 ), protostellar disks, etc. 
In the absence of a proof of the fluctuation dynamo at 
low Pm, the case for the dynamo origin of the small-scale 
magnetic fields in such systems (e.g., the observed fields 
in the solar photosphere) has been based on simulations 
done in the opposite regime of Pm > 1 [l(| ■ 

Previous numerical investigations of the onset of the 
fluctuation dynamo at low Pm ITt| revealed that 

the critical magnetic Reynolds number Rm c for dynamo 
increased with Re. In this Letter, we report that Rm c 
eventually reaches a finite limit as Re — > oo, i.e., we 
demonstrate for the first time that the fluctuation dynamo 
at asymptotically low Pm exists. The most important 
outstanding issue is whether the dynamo we have found 
is driven by the inertial-range motions at the resistive 
scale — if it is, its growth rate should be proportional to 
Rm 1 / 2 , which would make it a dominant (and universal!) 
field-amplification effect compared to any mean-field dy- 
namo due to the outer-scale motions. 
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FIG. 1: Growth/decay rate of (B 2 ) vs. Pm for 77 = 4 x 1CT 3 
(Rm ~ 60), ij = 2x 1CT 3 (Rm ~ 110), 77 = 10~ 3 (Rm ~ 230), 
77 = 5 x 10~ 4 (Rm ~ 450), and n = 2.5 x 10~ 4 (Rm ~ 830). 

Numerical Set Up. The equations of incompressible 
magnetohydrodynamics 

d t u + u • Vu = -Vp - u n \S/\ n u + B VB + f , (1) 
d t B + u- VB = B Vu + r;V 2 B (2) 

are solved in a periodic code of size 1, using a standard 
pseudospectral code. Here u is the velocity and B the 
magnetic field (in velocity units). The density is con- 
stant and equal to 1 . The velocity is forced by a random 
nonhelical homogeneous white-noise body force, which 
injects energy into the velocity components with wave 
numbers |k| < \/2ko, where ko = 2tt is the box wave 
number. Because the forcing is a white noise, the aver- 
age injected power is fixed: e — (u ■ f ) = 1, where the 
angle brackets stand for volume and time averaging. The 
numerical integration is continued only for as long as is 
needed to obtain a converged value of the growth/decay 
rate. In all of our runs, the magnetic field is energeti- 
cally much weaker than the velocity at all times, so the 
Lorentz force is never important. 

The maximum resolution we could afford was 512 3 . In 
order to increase the range of accessible Reynolds num- 
bers, we performed simulations both with the Laplacian 
viscosity (n = 2 in Eq. (p}) and with the 8th-order hy- 
perviscosity (n = 8). We define Rm = (u 2 ) 1 ^ 2 /nko, 
Rc = (u 2 ) 1 / 2 /vko, Pm = v/rj, where v — v% for the 
Laplacian runs and v = v c s — e/(|Vu| 2 ) for the hyper- 
viscous ones. We believe that using hyperviscosity is jus- 
tified for Pm <C 1 because the resistive scale in this limit 
is much larger than the viscous scale, 1^ 3> l v , and the 
magnetic properties of the system should be independent 
of the form of viscous regularization. 



Existence 0/ the Dynamo. In Fig. [T] we show the 
growth/decay rates of the magnetic energy (B 2 ) vs. Pm 
for five sequences of runs, each with a fixed value of 77. 
Thus, decreasing Pm is achieved by increasing Re while 
keeping Rm fixed. The growth rates arc calculated via a 
least-squares fit to the evolution of In ((-B 2 )) vs. time. We 
find that as Pm is decreased, the growth rate decreases, 
passes through a minimum and then saturates a constant 
value, i.e., at fixed Rm, 7(Rm,Re) — > 7oo(Rm) = const 
as Re — > 00. It is natural that such a limit exists because 
'77 >■ lv ~ LRe~ 3 ^ 4 — > and one cannot expect the mag- 
netic field to "know" exactly how small the viscous scale 
is. The nontrivial result is that, as Rm increases, the en- 
tire curve 7(Re, Rm) is lifted upwards, so the asymptotic 
values 7oo(Rm) > 0. Although we cannot at current res- 
olutions determine these positive asymptotic values, our 
judgement is that Fig. Q] provides sufficient evidence for 
claiming that such values exist and are positive. 

The Stability Curve. In Fig. [2j we show the recon- 
structed stability curve Rm c (Re) : each point on the curve 
is obtained by linear interpolation between a decaying 
and a growing case. We see that Rm c increases with 
Re, reaches a maximum Rm< max) ~ 350 at Re ~ 3000, 
and then decreases. Since the exact value of the vis- 
cous scale cannot matter when Re 3> Rm, Rm c (Re) — * 
Rm^°°' = const as Re — > 00. We cannot as yet obtain this 
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FIG. 2: The stability curve Rm c vs. Re. "Error bars" connect 
(Re, Rm) for decaying and growing runs used to obtain points 
on the stability curve. Stability curves based on the Laplacian 
and hyperviscous runs are shown separately. For comparison, 
we also plot the Rm c (Re) curve obtained in simulations em- 
ploying TGI [13 and TG2 forcing 111] (the three highest-Re 
points in the latter case were obtained by large-eddy simula- 
tions). The values of Re and Rm are recalculated according 
to our definitions, using the forcing wavenumber ko, rather 
than the dynamical integral scale as in [l^. [l7j. 



3 




FIG. 3: Normalized spectra of kinetic energy (compensated by fc 5//3 ) and (growing) magnetic energy for fixed n — 5 x 10 
(the Rm ~ 450 sequence from Fig. [1} and increasing Re: (a) Laplacian runs, (b) hyperviscous runs. 



asymptotic value. Unless the stability curve has multi- 
ple extrema (which we consider unlikely), Rm[°°' < 200. 
Note that this value is only about 3 times larger than 
the well-known critical value Rm c ~ 60 for the fluctua- 
tion dynamo at Pm > 1 0,0, [lj|. 

Considering that there is a small but measurable dif- 
ference between the stability curves for the Laplacian and 
hyperviscous simulations, how universal are our results? 
We believe it is plausible to argue, as we did above, that 
the existence of the dynamo in the limit Rm <§; Re — > oo 
does not depend on the nature of the viscous regulariza- 
tion. We also think that the asymptotic value Rm^ 00 ' 
is likely to be robust. However, the full functional de- 
pendence Rm c (Re) is certainly not universal. Indeed, let 
us consider what determines the shape of the stability 
curve in the transition region between the high- and low- 
Pm limits. When Re < Rm, l v > l n . As Re is increased, 
the spectral bottleneck associated with the viscous cutoff 
moves past the resistive scale l v until finally l v <C l v and 
the resistive scale is in the inertial range. This transition 
is illustrated by Fig.[3j The properties of the velocity field 
around the viscous scale are obviously not independent 
of the viscous regularization. Therefore, the 7(Re, Rm) 
and Rm c (Re) cannot be universal in the transition re- 
gion. In particular, since the bottleneck is narrower for 
the hyperviscous runs, the transition in the parameter 
space should be sharper. 

Magnetic- Energy Spectrum. The shape of the 
magnetic-energy spectrum is qualitatively different for 
Pm > 1 and Pm <C 1 (see Fig. [3]). At Pm above and 
just below unity, the spectrum has a positive slope and 
its peak is at the resistive scale. This is typical for the 
fluctuation dynamo at Pm > 1 — in the limit Pm 3> 1, 



the Kazantsev [HI fc +3 / 2 spectrum emerges 0j- As Pm 
is decreased, the spectrum flattens and then appears to 
develop a negative slope in the inertial range. At current 
resolutions, it is not possible to determine definitively 
what the asymptotic spectral slope is and whether the 
spectral peak is independent of Rm or moves with the 
resistive scale as fc Dca k oc Rm 3 / 4 . 

Simulations with a Mean Flow. 
2(il . 21 1 have been motivated 



Comparison with 
Several authors [13, [l 
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by the liquid-metal dynamo experiments to investigate 
the dynamo action at low Pm in numerical simulations 
where the forcing was spatially inhomogeneous and con- 
stant in time rather than random. The velocity field in 
these simulations consisted of a time-independent mean 
flow and an energetically a few times weaker fluctuating 
component (turbulence). The stability curves Rm c (Re) 
obtained in these studies have an entirely different origin 
than ours. In order to illustrate the difference, Fig. [2] 
shows the stability curves for simulations with Taylor- 
Green forcing, using published data 12|, |l7| . We see that 
the dynamo threshold for the simulations with a mean 
flow is much lower than for our homogeneous simulations. 
The difference is not merely quantitative. The mean 
flows in question are mean-field dynamos (even in the 
case of the nonhelical Taylor-Green forcing). This is con- 
firmed by the ordered box-scale structure of the growing 
magnetic field reported for these simulations at Pm > 1 
(the lower part of their stability curve). For Pm > 1, the 
threshold for the field amplification is Rm c ~ 10, which 
is a typical situation for mean-field dynamos [11]. The 
presence of magnetic energy at small scales is probably 
due to the random tangling of the mean field by turbu- 
lence, rather than to the fluctuation dynamo, because 
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Rm is below the fluctuation-dynamo threshold. The in- 
crease of Rm c with Re in these simulations is due to the 
interference by the turbulence with the dynamo proper- 
ties of the mean flow (l9l . |22| . It has not been checked 
whether the turbulence in these simulations might itself 
be a dynamo. Comparison of the two stability curves in 
Fig. [2] suggests that this can only happen at much larger 
Rm than studied so far. 

Outstanding Questions. The most important factual 
issue that remains unresolved is whether the low-Pm dy- 
namo we have found is due to the inertial-range mo- 
tions. The most important physical question is what is 
the physical mechanism that makes the dynamo possible. 

If the local (in scale space) interaction of the inertial- 
range motions with the magnetic field is capable of am- 
plifying the field, the dominant contribution to such a 
dynamo should be from the motions at the resistive scale 
l v — {jf/e) 1 ^ — LRm -3 / 4 [§], where the stretching 
rate is maximal. The growth rate of the magnetic en- 
ergy should scale as the stretching rate: 7 ~ (e/^) 1 ^ 2 ~ 
(Wrms/ XjRm 1 ' 2 . For Rm ^> 1, such a dynamo would 
always be faster than a mean-field or any other kind 
of dynamo associated with the outer-scale motions, be- 
cause the latter cannot amplify the field faster than at 
the rate ~ u Ims /L. Thus, the most pressing task for fu- 
ture numerical studies is to determine whether 7 scales 
as Rm ' or reaches an Rm- independent limit. If the lat- 
ter is the case, one will have to conclude that it is the 
outer-scale motions that act as a dynamo despite (or in 
concert with) the turbulence in the inertial range. Unlike 
the inertial-range dynamo, the characteristics of such a 
dynamo would not be universal. 

Theoretical predictions for a Gaussian white-noise ve- 
locity field (the Kazantsev [l8j model) strengthen the 
case for an inertial-range dynamo. For a certain range of 
scaling exponents of the velocity correlation function, it 
is possible to prove that the Kazantsev field is a dynamo 
23L 24 1 . The dynamo threshold for Pm 1 is predicted 
to be Rm[°°' ~ 400 (using our definition of Rm) — an 
overestimate by only (or by as much as) a factor of 2. The 
difficulty in deciding whether this theory applies lies in 
the unknown effect that assuming zero correlation time 
has on the dynamo properties of the inertial-range ve- 
locity field (in real turbulence, this correlation time is 
not only not small but also scale dependent). The main 
problem is that the result of Kazantsev theory is purely 
mathematical and that we do not have a physical model 
of the inertial-range dynamo. 

Conclusions. We have established that the fluctua- 
tion dynamo exists in nonhelical randomly forced homo- 
geneous turbulence of a conducting fluid with low Pm. 
The critical Rm in this regime is approximately 3 times 
larger than for Pm > 1. The nature of the dynamo and 
its stability curve Rm c (Re) are different from the low- 
Pm dynamo found in simulations and liquid-metal ex- 



periments with a mean flow. The physical mechanism 
that enables the sustained growth of magnetic fluctua- 
tions in the low-Pm regime is unknown. Is is not as yet 
possible to determine numerically whether the fluctua- 
tion dynamo is driven by the inertial-range motions at 
the resistive scale and consequently has a growth rate 
oc Rm ' , or rather is an outer-scale effect and has a 
constant growth rate comparable to the turnover rate of 
the outer-scale motions. 
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